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Abstract We review the observations of supernova rem- 
nants (SNRs) and pulsar-wind nebulae (PWNe) that 
give information on the strength and orientation of mag- 
netic fields. Radio polarimetry gives the degree of or- 
der of magnetic fields, and the orientation of the or- 
dered component. Many young shell supernova rem- 
nants show evidence for synchrotron X-ray emission. 
The spatial analysis of this emission suggests that mag- 
netic fields are amplified by one to two orders of magni- 
tude in strong shocks. Detection of several remnants in 
TeV gamma rays implies a lower limit on the magnetic- 
field strength (or a measurement, if the emission process 
is inverse- Compton upscattering of cosmic microwave 
background photons). Upper limits to GeV emission 
similarly provide lower limits on magnetic-field strengths. 
In the historical shell remnants, lower limits on B range 
from 25 to 1000 /iG. Two remnants show variability of 
synchrotron X-ray emission with a timescale of years. 
If this timescale is the electron-acceleration or radia- 
tive loss timescale, magnetic fields of order 1 mG are 
also implied. In pulsar-wind nebulae, equipartition ar- 
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guments and dynamical modeling can be used to in- 
fer magnetic-field strengths anywhere from ~ 5 /zG 
to 1 mG. Polarized fractions are considerably higher 
than in SNRs, ranging to 50 or 60% in some cases; 
magnetic-field geometries often suggest a toroidal struc- 
ture around the pulsar, but this is not universal. Viewing- 
angle effects undoubtedly play a role. MHD models of 
radio emission in shell SNRs show that different orienta- 
tions of upstream magnetic field, and different assump- 
tions about electron acceleration, predict different ra- 
dio morphology. In the remnant of SN 1006, such com- 
parisons imply a magnetic-field orientation connecting 
the bright limbs, with a non-negligible gradient of its 
strength across the remnant. 



1 Shell supernova remnants: review 
1.1 Introduction 

Supernova remnants and pulsar- wind nebulae are promi- 
nent Galactic synchrotron sources at radio and, often, 
X-ray wavelengths. The spatial and spectral analysis of 
the synchrotron emission can be used to deduce or con- 
strain magnetic-field strengths and orientations. In this 
review, "supernova remnant" (SNR) will be used for 
shell (i.e., non-pulsar-driven) remnants; if the emission 
is due to a pulsar, whether the object is young (like 
the Crab Nebula) or old (as in Ha bow-shock nebulae) 
the object will be termed a pulsar- wind nebula (PWN). 
Thus a SNR may contain a PWN, a combination some- 
times called a "composite" SNR. 

SNRs are primarily radio objects. In the Milky Way, 
274 are listed in Green's (2009) well-known catalogue 
(http://www.mrao.cam.ac.uk/surveys/snrs/); most are 
known only by radio emission. SNRs are well-studied 
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only in Local Group galaxies where they are bright and 
large enough to be well-imaged with radio interferome- 
ters. Galactic remnants range in angular size from less 
than 2' to many degrees, and represent a range in ages 
from about 100 years (Gl.9+0.3) to over 10 5 yr, though 
age estimates are difficult for very old objects. 

1.2 SNR Dynamics 

SNRs evolve through various phases as their shock waves 
decelerate in the ambient medium. Very soon after an 
explosion (either Type la or core-collapse) deposits ~ 
10 51 erg into the ISM, a shock wave is driven into 
the surrounding circumstellar medium (CSM: modified 
by the progenitor) or unmodified interstellar medium 
(ISM). Core-collapse SNe eject several solar masses with 
a range of velocities, typically of order 5,000 km s — 1 . 
Type la's are thought to represent the thermonuclear 
disruption of a white dwarf, so eject 1.4 M© at about 
10,000 km s _1 . Not long after the explosion, decelera- 
tion of the outer blast wave causes the formation of a 
reverse shock as inner ejecta are forced to decelerate. 
This begins the ejecta phase of evolution, characterized 
by the presence of both the forward shock (blast wave) 
heating ISM or CSM, and the reverse shock heating 
ejecta. Cooling timescales are much longer than dynam- 
ical timescales, so this phase is essentially adiabatic. 
Depending on the density structure in the ejecta and 
surroundi ng material, th e shock radius may evolve as 
t°- 6 -t 9 t)Chevalierlll982l ). Eventually the reverse shock 
moves in to the center of the SNR and, after what may 
be extensive reverberations, disappears. By this time 
(several thousand years for typical parameters), sev- 
eral times the ejected mass have been swept up by the 
blast wave, and the remnant is settling into the Sedov 
phase, well described by the similarity solution for a 
point explosion in a uniform medium, with shock ra- 
dius obeying R oc t° . This phase is adiabatic as well. 
Eventually, cooling times in shocked ISM become com- 
parable to the remnant age, and shocks become radia- 
tive: lossy and optically prominent. This typically oc- 
curs when shock velocities v s drop to around 200 km 
s . Denser regions around the remnant periphery may 
become radiative before other parts; even young SNRs 
(like Kepler, SN 1604) can show radiative shocks in 
some regions. Once the bulk of the blast wave is ra- 
diative, the deceleration is more marked. The remnant 
interior will remain hot for some time, producing what 
is known as a "pressure-driven snowplow." Eventually, 
the remnant becomes sufficiently confused as to lose its 
identity, and the remaining kinetic energy is dissipated 
as sound waves i n the ISM. Various b l ast-wa ve solutions 
are described in Ostriker fc McKee! (|l988h : particular 



applic ations to SNRs are made in Truelove fc McKeel 
(19991). 



1.3 Radio inferences 

The prevalence of radio emission from adiabatic-phase 
SNRs, with low shock compression ratios, indicates that 
particle acceleration must be ongoing; simple compres- 
sion of ambient Galactic magnetic field and cosmic-ray 
electrons would produce both in sufficiently brigh t emis- 
sion and the wrong spectrum ( Reynolds 20081) . Since 
1977, the traditional explanation for this electron ac- 
celeration has been diffusive shock acceleration (DSA; 
see Blandford & Eichler 1987 for a comprehensive re- 
view). Electrons emitting synchrotron radiation at ra- 
dio wavelengths have energies of order 1-10 GeV 
(E = U.7y/iy(GRz)/B(fiG) GeV). The observed ra- 
dio power-law spectra S v oc v~ a ^ with a ~ 0.5 (see 
Figure Q]), imply from synchrotron theory a power-law 
energy distribution of electrons KE~ S electrons cm' 3 
erg -1 with s = 2a + 1 ~ 2. Radio observations do 
not allow the separate deduction of energy in electrons 
and in magnetic field; the synchrotron emissivity of this 
power-law distribution of electrons is proportional to 
_£Cf?( s + 1 )/ 2 , or roughly to the product of energy den- 
sity in electrons and in magnetic field. If particle and 
magnetic-field energies were in equipartition (the minimum- 
energy state producing an observed spectrum), the magnetic- 
field strength could be deduced. For a spherical remnant 
of radius R pc at distance D kpc with a = 0.5, 

B = 20 ((1 + k)SvDl pc /{4>Rl c )f 7 f,G, (1) 

where Sg is the flux density at 1 GHz in Jy, 4> is the 
volume filling factor, and k is the ratio of energy den- 
sity in ions to that in electrons (e.g., Pacholczyk 1970). 
But there is no obvious reason that equipartition should 
hold: SNRs are typically very inefficient synchrotron ra- 
diators, with a very small fraction of the 10 51 erg re- 
quired to produce the observed emission. Thermal en- 
ergies are very much larger than either particle or field 
energy. Thus equipartition magnetic-field strengths are 
really just a measure of (approximately) remnant sur- 
face brightness. 

The question of magnetic-field amplification in strong 
SNR shocks, to be discussed at length below, is inti- 
mately related to the question of shock modification 
by cosmic-ray pressure (efficient shock acceleration, or 
nonlinear DSA). Thus evidence for nonlinear DSA is 
related to estimation of magnetic-field strengths. One 
prediction of this theory is the gradual deceleration of 
incoming flow (in the shock frame) by forward-diffusing 
cosmic rays, so that the compression ratio r varies con- 
tinuously from a possibly large value (far upstream ve- 
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SNR Spectral Indices 
(Green 2009) 



dynamical precursor 




Fig. 1 Histogram of shell SNRs with fairly well-measured radio 
spectral indices, from Green 2009. PWNe are excluded. 



locity/downstream velocity, perhaps 10 or more) to about 
3, at which point most calculations show the develop- 
ment of a "thermal subshock" on a length scale of ther- 
mal proton gyroradius, much smaller than other length 
scales. (See Figured) In test-particle (inefficient) DSA, 
well-known results give the spectrum of accelerated par- 
ticles as a power-law with energy index s depending 
only on the shock compression ratio, s = (r + 2)/(r — 1). 
Since strong adiabatic shocks are expected to have r = 
4 (for adiabatic index 7 = 5/3), we predict s — 2 which 
implies a = 0.5, close to the mean of the distribution 
of over 200 SNR values from Green's catalogue. How- 
ever, the dispersion in that distribution is substantial 
(see Figure [1]); explaining values of a of 0.6 - 0.7, typi- 
cal for young SNRs, or as low as 0.3, sometimes seen in 
older SNRs, is not straightforward. Low Mach-number 
shocks will have lower compressions: in fact, for Mach 
number M, r" 1 = (7- l)/( 7 + 1) + 2/(( 7 + l)M 2 ). For 
A4 < 10, this can significantly increase the steepness 
of the predicted spectrum; but this is demonstrably 
not the case in the young remnants with the steeper 
spectra. No variation in shock Mach number can ex- 
plain a < 0.5; traditional explanations involve confus- 
ing thermal emission or other instrumental effects, but 
remain unconvincing. 

A consistent explanation for steeper radio spectra 
in young remnants is available, however. Particles of 
a particular energy will diffuse (on average) a certain 
distance ahead of an efficient shock, where they will 
see a particular effective compression ratio, which fixes 
the slope of the distribution function near that energy. 
If, as expected, more energetic particles diffuse further, 
they will see larger compression ratios and the spec- 
trum will become flatter. Thus one predicts spectra 




Larger 

compression 

ratio 



upstream 



1=0 
shock 



downstream 



Fig. 2 Schematic velocity profile of a shock wave with upstream 
velocity u\ and downstream U2- Dotted line: test-particle shock 
(velocity discontinuity). Solid line: shock modified by cosmic rays. 



flattening (hardening) to higher energy, an effect first 
pointed out by Eichler (1979). This effect was quanti- 
fied in a Monte Carlo simulation of electron accelera- 
tion, and compared with radio observations of young 
SNRs, which did seem to show the effect (Ellison & 
Reynolds 1991; Reynolds & Ellison 1992; see Figure 3). 
The prediction for the emitted synchrotron spectrum 
depends on the magnetic-field strength; stronger fields 
mean lower-energy particles producing radiation in a 
particular radio-astronomical bandpass. Reynolds & El- 
lison found that explaining the relatively steep spectra 
of Tycho's and Kepler's SNRs required magnetic fields 
of up to 1 mG and above, though with substantial er- 
rors. 

The degree of order in the magnetic field in SNRs 
can be directly inferred from the degree of polariza- 
tion of synchrotron radiation (subject to observational 
issues such as Faraday depolarization and resolution ef- 
fects [beam depolarization] ) . The degree of linear polar- 
ization of power-law synchrotron radiation from elec- 
trons in a uniform magnetic field is given by P = (s + 
l)/(s + 7/3), so that s = 2 =$> P = 69%. Such high val- 
ues are never observed in SNRs; a few older SNRs such 
as DA 530 show P 



50% in some regions (jLandecker et al 
1999, Figure 0]) , but historical shells show much lower 
values, 10 - 15% (see references in Reynolds & Gilmore 
1993). Very roughly, the observed polarized fraction 
over the maximum possible gives an estimate of the 
fractional energy in th e ordered component of the mag- 
netic field ( Burn 1966[) ; thus in Tycho, for instance (see 
Figure [5] [Reynoso et al. 1997]), less than 15% of the 
magnetic energy is in an ordered component. However, 
the direction of that component is radial, a peculiar 
property shared by the other historical shells SN 1006, 
Cas A, Kepler, and Gil. 2-0. 3. The nature of this ra- 
dial component is unclear; while radial motions such as 
predicted to occur at the Rayleigh- Taylor unstable con- 
tact discontinuity between shocked ejecta and shocked 
ISM might enhance radial field, this effect should occur 
somewhat interior, while radial fields are observed im- 
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Fig. 3 Observed total flux observations of Tycho's and Kepler's 
SNRs and SN 1006, with model spectra superposed. All three 
spectra derive from the same electron spectrum but with var- 
ied magnetic field. For Tycho and Kepler, the dashed line corre- 
sponds to B = 10 mG, the solid line to 1 mG, and the dot-dashed 
line to 0.1 mG; for SN 1006, the dashed, solid, and dot-dashed 
lines correspond to B = 10 mG, 0.1 mG, and 1 fiG, respectively. 
The dotted line shows the spectrum that would be produced if 
electrons had the same energy distribution as protons, showing 
how differently the nonlinear shock treats electrons and protons. 
(Reynolds & Ellison 1992) 



mediately at remnant edges (e.g., Tycho: Dickel et al. 
1991). Various MHD simulations have been performed 
in an attempt to unders tand this effect (jJun fc Norman 
1996tlJun fc Joneslll999h . In older remnants such as DA 
530, magnetic fields are generally confused or tangen- 
tial. Tangential fields can readily be explained by simple 
compression in radiative shocks with large r, and this 
is the conventional explanation. 



1.4 X-ray and gamma-ray inferences 

The discovery of X-ray synchrotron emission from SNRs 
dates from the early 1980's (Toor 1980; Becker et al. 1980; 
Reynolds & Chevalier 1981) but observational and the- 
oretical confusion prevented its wide acceptance un- 
til the unmistakable observations with ASCA of the 
remnant of SN 1006, demonstrating conclusively that 
central emission was thermal while the bright limbs 
showed lineless spectra well described by power laws 
(Koyama et al. 1996). (See Reynolds 2008 for a de- 
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Fig. 4 Radio contours and polarized B-vectors at 1420 MHz for 
DA 530 (Landecker et al. 1999). Direction of the sky-plane com- 
ponent of the magnetic field is at right angles to the directions of 
the i?-vectors, i.e., predominantly tangential in the bright limbs. 
Foreground Faraday rotation has been corrected for, though it is 
small. 
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RIGHT ASCENSION (J2000) 

Fig. 5 Polarization vectors in Tycho's SNR at 1375 MHz 
(Rcynoso et al. 1997). These have been corrected for foreground 
Faraday rotation and rotated by 90° to show the direction of the 
magnetic field. 



tailed review of high-energy nonthermal emission from 
SNRs.) At this time, four Galactic SNRs show X-ray 
spectra dominated by synchrotron emission: Gl.9+0.3, 
SN 1006, G347.3-0.5 (RX J1713. 7-3946), and G266.2- 
1.2 ("Vela Jr." or RX J0852-4622). Of these, the first 
two are young, with symmetric X-ray morphologies, 
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and show thermal X-ray emission from fainter regions. 
SN 1006 is widely accepted to be the remnant of a Type 
la supernova (evidence includes the presence of Balmer- 
dominated optical emission from the forward shock, and 
the high Galactic latitude of 15°); the evidence is not as 
firm for Gl. 9+0.3, but it is also su spected to be a Type 
la remnant ( Reynolds et al. 20081) . The other two have 
much larger angula r sizes and considerably more irreg- 
ular morphologies ( Slane et al. 19991 2001 ): they also 
show no trace of thermal X-rays, putting severe upper 



limits on the presence of any thermal gas ([Ellison et al 



2010l) . Both also contain X-ray point sources, suggest- 
ing a core-collapse origin. In addition to these four ob- 
jects, all other historical shells (Tycho = SN 1572, Ke- 
pler = SN 1604, Cas A (unseen SN around 1680), G11.2- 
0.3=SN 386, and RCW 86 = SN 185?) show synchrotron 
emission in regions, typically (but not always) "thin 
rims" at the remnant periphery. All, in addition to 
SN 1006, also show hard continua in the integrated 
X-ray flux seen by the (non-im aging) RXTE satellite 
(jAllen. Gotthelf. fe Petrell999h . Detections above 8 keV 
have also been reported by INTEGRAL (Cas A, Re- 
naud et al. 2006; SN 1006, Kalemci et al. 2006; G347.3- 
0.5, Krivonos et al. 2007) and by non-imaging instru- 
ments BeppoSAX (Cas A to 80 keV; Favata et al. 1997) 
and Suzaku HXD (Cas A [to 40 keV], Maeda et al. 2009; 
G347.3-0.5 [to 40 keV], Tanaka et al. 2008; and Tycho 
[to 30 keV], Tamagawa et al. 2009). Arguments that 
the radiation is synchrotron essentially turn on the ex- 
clusion of all other possible processes (see arguments 
summarized in Reynolds 2008). Synchrotron photons 
at 4 keV energy (v ~ 10 18 Hz) imply electron energies 
of order 100(^/10 fiG) 1 / 2 TeV. 

Power-law fits to the X-ray spectra are much steeper 
than radio spectra, indicating that the X-rays come 
from the cutting-off tail of the electron distribution. 
Three limitations might be imagined for the maximum 
energy to which electrons can be accelerated due to 
DSA: finite age or size of the remnant, lack of MHD 
scattering waves above some A(max), or radiative losses. 
The first two mechanisms would restrict ion accelera- 
tion as well, but radiative losses would affect electrons 
only. The cutoffs are likely to be exponentials or mod- 
ified exponentials, so that the electron distribution is 
roughly given by N(E) — KE~ s e~ E / Era ^. In a uni- 
form magnetic field, such a distribution would give rise 
to a spectrum S v oc v~ a e~^j v l Vc approximately, drop- 
ping off rather slowly. (A more careful calculation in the 
case of synchrotron losses predicts a spectrum with the 
same exponential factor e~ v 7 "/^ ; Zirakashvili & Aharo- 
nian 2007). Here v c cx E 2 lax B is the "rolloff" frequency 
corresponding to the peak frequency emitted by elec- 
trons with energy £ mM . 



In the cutoff part of the spectrum, we observe the 
competition between acceleration and loss rates. In stan- 
dard DSA, the time for a particle to reach a relativistic 
energy E (from E.- L <C E) is given by t(E) ~ k/u 2 . 
where k = A m f p c/3 is the diffusion coefficient. It is of- 
ten assumed, but not demanded by the data, that k is 
"Bohm-like," that is, that the mean free path A m f p is 
proportional to the particle gyroradius r g = E/eB (for 
extreme- relativistic energies). If we write A m f p = nr g , 
with rj the "gyrofactor," depending on the amplitude of 
scattering turbulence, then r\ = 1 is Bohm diffusion and 
r\ =const.> 1 is "Bohm-like" . In commonly used "quasi- 
linear theory," rj(E) = (SB/B)~ 2 where SB is the am- 
plitude of MHD turbulence resonant with particles of 
energy E. (So r\ = const, means a particular spectrum 
of turbulence, equal energy per unit logarithmic band- 
width). Bohm diffusion is often taken as a limiting case, 
SB ~ B, but it has been argued that SB ^> B is the 
possible outcome of some cosmic-ray driven instabili- 
ties (Bell & Lucek 2001; Bell 2004). Relevant for using 
X-ray spectra as a diagnostic of magnetic fields is that 
high fields give short mean free paths, small diffusion 
coefficients, rapid acceleration, and higher maximum 
particle energies. 

Two classes of argument for high magnetic fields in 
young SNRs come from the last decade of observations 
with high-resolution, high throughput X-ray observato- 
ries: Chandra, XMM-Newton, and Suzaku. A morpho- 
logical argument is based on the commonly seen phe- 
nomenon of "thin rims," in which synchrotron X-ray 
emission occurs at remnant peripheries in very narrow 
tangential features coincident with the shock location 
as inferred, e . g., fro m H a observations (see Fi gure Et. 
Bamba et al. ( 2003| ) and Vink fe Laming ( 2003 ) argued 
that the small radial extent resulted from synchrotron 
losses on electrons, to infer strong amplification of B: 

B ~ 200(u s /1000 km s _1 ) 2/a (ij;/0.01 pc)~ 2/3 G (2) 



( Parizot et al.l 120061 ) . where w is the filament width. 
Values of B from 60 to 230 fiG have been inferred 
for Tycho, Kepler, SN 1006, Cas A, and G347.3-0.5 



( Parizot et al. 20061) . These field strengths apply only 



to regions where X-ray synchrotron emission is seen, of 
course. An alternative explanation for the disappear- 
ance of synchrotron emission such a short distance be- 
hind the shock is that the magnetic field i s primarily 



turbu lent, and decays on a short lengthscale (jPohl. Yan. fe Lazarian 
20051) . Tests attempting to discriminate between these 
two possibilities somewhat favor th e radiative-loss ex- 



plana tion, but are not yet conclusive (jCassam-Chenai et al 
2007) . One particular difficulty is the presence of thin 



rims in radio images as well (see Figure [7]); since radio- 
emitting electrons have enormously longer loss timescales, 
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Fig. 6 Raw Chandra counts image of the NE region of the limb 
of SN 1006, showing "thin rims" of synchrotron radiation (Long 
et al. 2003). The brightness contrast of the very sharp rims is 
about a factor of 2. The inner rim has a a clear He* counterpart 
[Winkler et al. 2003]. 
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Fig. 7 Left: Chandr a image of Tycho's SN R (CXC). Right: VLA 
linage at 1420 MHz llRevnoso et al.ll 199^ . Note the thin rims in 
the X-ray image (e.g., NE quadrant), with some thin radio rims 
as well. 



disappearance of the magnetic field is the only alter- 
native. It is possible that both processes operate; theo- 
retical models for strong magnetic-field amplification in 
modified shocks produce turbulent magnetic field which 
might well have some decay mechanism to compete with 
electron energy losses. 

Another argument for (locally) strong magnetic fields 
is the discovery of variations in brightness on timescales 
of a few years of small features in X-ray synchrotron 
emission m Cas A (|Patnaude fc Fesenl2007l) and G347.3- 



0.5 (jUchivama et al.l 120071 ). Again, factors of 2 drop in 
brightness on a timescale of years require B ~ 1 mG. 
Brightness increases on similar timescales are also ob- 
served; demanding acceleration times of years gives sim- 
ilar estimates for B. However, very strong levels of mag- 
netic turbul ence also naturally predict these kinds of 
fluctuations ( Bvkov. Uvarov. fc Ellisonl2008 ; Bvkov et al 



2009) without necessarily requiring quite as high values 
of B lnls . However, the remnant of SN 1006 shows no 



such small features or brightness variations (jKatsuda et al 
2010), perhaps as a result of a Type la origin with 
more-or-less uniform surrounding material. (The ab- 
sence of short-term variations does not imply absence of 
magnetic-field amplification; very high B values could 
be generated in a quasi-steady state behind the shock, 
without producing short-term brightness fluctuations.) 

New arguments relative to magnetic-field strengths 
in SNRs have resulted from the new generation of GeV- 
TeV observational capabilities. The Fermi Gamma-Ray 
Space Telescope has been mapping the sky between 
about 0.2 and 300 GeV since its launch in 2008, with 
gradually improving statistics on any given region of 
sky. One important early re sult is the detectio n of Cas 
A between 0.5 and 50 GeV (|Abdo et al.ll2010al see Fig- 
ure [5]). The emission may be either leptonic, that is, 
electron bremsstrahlung plus inverse-Compton upscat- 
tering of local photon fields, or hadronic, from the decay 
into photons of ir° mesons produced in inelastic colli- 
sions between cosmic-ray protons and local thermal gas. 
The bremsstrahlung and hadronic contributions require 
knowledge of the gas density; the IC contribution re- 
quires knowledge of appropriate photon fields. The lep- 
tonic model shown in Figure [5] assumes a mean density 
of 26 cm -1 , which requires a mean magnetic field (aver- 
aged over the entire emitting region of Cas A) of about 
0.12 mG. If the actual emission is hadronic, the IC con- 
tribution must be lower (fewer relativistic electrons), 
demanding a larger magnetic field to produce the ob- 
served synchrotron fluxes at radio energies. (This argu- 
ment was originally made based on early upper limits 
above 100 MeV from COS-B and SAS-2; Cowsik and 
Sarkar [1980] derived a lower limit of about 1 mG.) 

Similarly, any TeV detection or upper limit places 
lower limits on the mean magnetic field in regions con- 
taining relativistic electrons. So far, four SNRs have 
been imaged with the HESS air-Cerenkov telescope ar- 
ray in Namibia. These include three of the synchrotron- 
dominated remnants: SN 1006, G347.3-0.5, and Vela Jr. 
(The fourth, Gl.9+0.3, is too small for HESS imaging 
and too faint for a constraining upper limit). In addi- 
tion, the (prob ably) historical shell S NR RCW 86 has 
been detected ( Aharonian et al. 20091 ). In all cases, the 
TeV morphology tracks the X-ray morphology surpris- 
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Fig. 8 Spectrum o f GeV emission from Cas A as seen by Fermi 
llAbdo et aLlbOlOaT) . The solid curves are leptonic models (dots, 
IC; dashes, bremsstrahlung) with magnetic-field values shown. 
If the emission is actually hadronic, the magnetic-field strengths 
must be larger to suppress the IC component. 



ingly closely. However, detailed attempts to model the 
emission as either due to hadronic or to leptonic pro- 
cesses run into significant problems. For hadronic mod- 
els, required target gas densities seem to exceed mea- 
surements or upper limits on local thermal gas, while 
simple leptonic models do not do a good job reproduc- 
ing the spectral shape and may also require implausibly 
small filling factors of magnetic field. However, lower 
limits to (B) do n ot depend much on the detailed mod- 
eling. In SN 1006 (|Acero et al.ll2010l ). B > 30 tiG is re 



quire d; for G347.3-0.5, a leptonic model (jLazendic et al 



20041 ) demands a filling factor of magnetic field of only 
1% to avoid overpredicting synchrotron emission. Within 
that 1% of the volume, the field is a modest 15 /xG. 
(These estimates are relatively insensitive to model de- 
tails; even a model that does not use DSA at all as 
a particle acceleration mechanism, but rather stochas- 
tic acceleration in the downstream region, fits the most 
recent GeV and TeV data for G347.3-0.5 with a lep- 
tonic model with B ~ 12 riG; Fan, Liu, & Fryer 2010). 
Hadronic models typically invoke much higher fields 
(Zirakashvili & Aharonian [2010] find B = 127 fiG in 
G347.3-0.5). 

HESS has also detected several other sources coin- 
cident with shell SNRs, such as CTB37B (HESS J1713- 
381; Aharonian et al. 2008) and HESS J1731-347 (Acero 
et al. 2009). These may be further members of the 
nonthcrmal-X-ray class, older SNRs interacting with 
molecular clouds, or PWNe. They await further obser- 
vational clarification. 

The question of magnetic-field amplification is inti- 
mately connected with that of efficient particle accel- 
eration through the proposals of Bell & Lucek (2001) 
and Bell (2004) that cosmic-ray driven instabilities can 



Fig. 9 SN 1006 with Chandra (NASA/CXC). Red, 0.50 - 0.91 
keV; cyan, 0.91 - 1.34 keV; blue, 1.34 - 3.00 keV. Blue-white rims 
are line-free synchrotron emission; red interior is dominated by 
oxygen ejecta. 



greatly increase the magnetic-field strength. Evidence 
for high magnetic fields is then taken as indirect evi- 
dence for an energetically significant component of cos- 
mic rays (which must of necessity be ions) accelerated in 
the shock. This argument can be reversed: if evidence is 
found for efficient ion acceleration, then magnetic fields 
are likely to be amplified. 



1.5 Obliquity dependence 

The nature of magnetic-field amplification at strong 
shocks is a subject of intense theoretical activity. Ob- 
servations can provide powerful constraints on some as- 
pects of such amplification. In particular, two SNRs, 
SN 1006 and Gl.9+0.3, show very simple bilaterally 
symmetric morphology which may contain clues to the 
dependence of particle acceleration and magnetic-field 
amplification processes on the shock obliquity angle #Bn 
between the mean upstream magnetic field direction 
and the shock velocity. Section [3] discusses in detail a 
model for the radio morphology of SN 1006, as rep- 
resentative of this bilaterally symmetric class of SNRs 
(BSNRs). 

The most natural explanation for the bilateral sym- 
metry evident in Figures [9] and [TT] is the variation in #Bn 
one would encounter at different locations as a spheri- 
cal shock encounters a roughly uniform magnetic field. 
If the field lies in the plane of the sky, the range in 
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Fig. 10 RollofE frequencies measured in 10" deep (radially; much 
larger in the azimuthal direction) rectangular regions around the 
outer edge of SN 1006 (red points: an interior rim) (Katsuda et 
al., in preparation). The exposure time is much longer in the NE 
region, so boxes are smaller. 



Fig. 12 Variations of u TO n with radius and azimuth, for 
Gl.9+0.3. Red: "Ears" extending beyond the main shell to E 
and W. Blue: main shell. Green: Just interior to bright regions of 
main shell. Note the similar trend at much poorer resolution to 
that seen in Figure [TO] 




Fig. 11 Gl.9+0.3, the Galaxy's youngest SNR (Chandra; 
Borkowski et al. 2010). Red, 1-3 keV; green, 3-5 keV; blue, 
5 — 8 keV. Except for faint thermal emission in the center and 
north, the emission is dominantly synchrotron. 



A simple model adequately characterizing most syn- 
chrotron X-ray emission from shell SNRs is the "sr- 
cut" model in the XSPEC software package, which is 
the spectrum emitted by a power-law electron distribu- 
tion with an exponential cutoff. The frequency corre- 
sponding to the cutoff energy i5 max is v r o\\, the "rolloff" 
frequency (so called because the synchrotron spectrum 
from this distribution actually steepens rather gradu- 
ally). v TO \\ serves as a good indicator of the maximum 
energy reached by electrons at different locations, al- 
though depending on the mechanism limiting that max- 
imum energy, v m \\ has different dependencies on physi- 
cal parameters. 

If the acceleration rate of particles varies with shock 
obliquity (independently of any obliquity-variations of 
the magnetic field), we can parameterize such effects in 
terms of a quantity Rj: i?j(6* Bn , 77, r) = T(6> B n)/T(#Bn = 
0). (This means that the value of -q is that for O-Qn = 
0; any obliquity-dependence is subsumed into Rj.) We 
scale to typical values for young SNRs: 1*3.5 = u s h/3000 
km s _1 ; t 3 = </1000 yr; B 10 = 5/10 /iG; and A17 = 
Amax/10 17 cm. T hen the values of v rr A\ for each limiting 
mechanism obey (|Revnoldsll2008l) 



obliquities would be to 7r/2; if there is a line-of-sight 
component, the range drops until it is zero for B along 
the line of sight. Later in this chapter we describe a de- 
tailed calculation of the radio synchrotron morphologies 
resulting from different assumptions about the nature 
of the upstream magnetic field and of particle acceler- 
ation. Here we focus on the X-ray spectral variations 
each remnant displays. 



hv m \\{&ge) 
hv ro \\{esc) 

hv ro \\(].OSS) 



I 4 t 2 



(rjRj)- 2 keV 



2Bf 



2 keV 



2«8,Mi)"' keV. 



(3) 
(4) 
(5) 



The effective limiting mechanism is the one producing 
the lowest value of -E max ; this may vary around the 
remnant periphery. Note that if radiative losses limit 
electron acceleration, the resulting v m \\ is independent 
of the magnetic-field strength. 
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Figure [TU] shows the variation in u m \\ around the 
periphery of SN 1006 (Katsuda et al., in preparation); 
similar data for Gl.9+0.3 (Fi g. 11) at much lower res- 
olution ( Reynolds et al J [20091) are shown in Figure [T"2l 
Variations by large factors occur. Equations 3-5 above 
show that achieving variation by over an order of mag- 
nitude is not easy, especially for loss-limited acceler- 
ation. Even in the absence of any cosmic-ray-driven 
magnetic-field amplification, one would expect a vari- 
ation in B by about the shock compression ratio (4 
or more) simply due to compression of the tangential 
component of B due to flux-freezing. While there are 
no comprehensive predictions, calculations of magnetic- 
field amplification normally assume parallel shocks (#Bn ~ 
0), with the process being less effective at higher obliq- 
uities. The strong independence for age or escape- limited 
acceleration could easily explain the observations - but 
with the serious side effect of predicting that the result- 
ing maximum energies (-E max ~ 10 - 100 TeV) would 
apply to ions as well as electrons, putting at least these 
two SNRs out of the running for producing cosmic rays 
up to the "knee" , the slight steepening in the integrated 
spectrum of Galactic cosmic rays at Earth. Loss-limited 
maximum energies do depend on the shock speed, but 
the required variation by over 30 in v TO \\ would demand 
an implausible factor of 5 varia tion in the cu rrent shock 
velocity. Anisotropic diffusion (|jokipiilll9871 ) could pro- 
duce more rapid acceleration at quasi-perpendicular shocks 
than quasi-parallel (Rj{0Bn) < Rj(0)), easily explain- 
ing the required amount of variation (Figure II 3|) . How- 
ever, most workers believe that this effect may not oc- 
cur in highly turbulent media. In addition, electron in- 
jection into the DSA process may be more difficult at 
quasi-perpendicular shocks. As we show later in this 
chapter, the radio morphology may be more accurately 
reproduced with a quasi-parallel model, leaving the ex- 
planations of Figures [TUJ and [T^] as an open question. 



1.6 Summary 

We summarize here the inferences on magnetic-field 
strength and geometry described above. 

1. From radio observations, equipartition values of mag- 
netic field strength are in the ~ 10 /iG range, but 
there is little physical motivation to assume equipar- 
tition. 

2. Radio polarization studies show that in young SNRs, 
the magnetic field is largely disordered, with a small 
radial preponderance. In older, larger SNRs, the field 
is often disordered but sometimes tangential. 

3. Curvature (spectral hardening to higher frequency) 
is observed in the radio spectra of Tycho and Kepler. 
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compression ratio r = 4 
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Fig. 13 Variation of v IO \\ with azimuth (Reynolds 1998), ac- 
cording to the anisotropic diffusion model of Jokipii (1987). The 
required amount of variation of u vo \\ can be produced, but only 
if the shock is perpendicular at the bright limbs. B\ and B2 are 
the upstream and downstream magnetic-field strengths, and G is 
the ratio of an electron's energy-loss time in a field B\ to that in 
the actual obliquity-dependent combination of B\ and B2 as it 
scatters back and forth across the shock; see Reynolds 1998. 



A nonlinear shock acceleration model can explain 
this with magnetic field strengths of 0.1 - 1 mG 
(average over the emitting regions). 
Thin rims of X-ray synchrotron emission in a few 
young remnants require B ~ 50 — 200 /iG in the 
rims, if they are due to synchrotron losses on down- 
stream-convecting electrons. However, thin radio rims 
are sometimes seen as well; they require that the 
magnetic field disappear somehow, presumably be- 
cause it is a wave field which damps. 
Brightening and fading of small X-ray synchrotron 
features in G347.3-0.5 and Cas A require B ~ 1 mG, 
if they represent acceleration and loss times for elec- 
trons. Fields smaller by a factor of several are pos- 
sible if the fluctuations are due to strong magnetic 
turbulence. 

Large azimuthal variations in the rolloff frequency 
in SN 1006 and Gl.9+0.3 are difficult to explain for 
a conventional picture of loss-limited acceleration in 
parallel shocks. 

For Cas A, the detection at GeV energies with Fermi 
requires B > 0.1 mG to avoid overproducing the 
GeV emission with electron bremsstrahlung. 
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8. TeV emission seen in four shell SNRs is not well 
explained by either leptonic or hadronic processes. 
However, if it is hadronic, the magnetic fields im- 
plied are of order 100 /zG, while leptonic models 
require much lower fields. 



2 Magnetic fields in pulsar-wind nebulae 

2.1 Introduction 

All young pulsars are slowing down in their rotation. 
The corresponding rate of change of rotational kinetic 
energy is enormous, typically in the range E = 10 32 — 
10 39 ergs s" 1 . In most cases only a negligible amount 
of energy is contained in the electromagnetic radia- 
tion corresponding to the neutron star's observed pulsa- 
tions. Rather, the bulk of the rotational power is usually 
deposited into a relativistic magnetized particle wind 
that flows outward from the pulsar. 

External pressure will cause this wind to abruptly 
decelerate at a termination shock. Beyond the termina- 
tion shock, the pulsar wind thermalizes in pitch angle 
and radiates synchrotron emission, resulting in a pulsar 
wind nebula (PWN). The best known PWN is the Crab 
Nebula (Fig. HP!]) , powered by the central young pul- 
sar B0531+21. The Crab Nebula radiates synchrotron 
emission across the electromagnetic spectrum. In X- 
rays, its extent is small, reflecting the relatively short 
synchrotron lifetimes of electrons emitting at these en- 
ergies. The optical synchrotron emission is larger, re- 
flecting the longer lifetimes of the corresponding elec- 
tron population. And finally, the radio emission shows 
the full extent of the source, with a radiative lifetime 
exceeding the age of the source (~ 1000 years). 



2.2 Observational Properties 

PWNe and SNRs have a variety of common proper- 
ties: they both result from a supernova explosion, and 
they both are characterized by their polarized radio 
synchrotron emission. This has led to considerable con- 
fusion in the literature. For example, the Crab Nebula 
is often described as a "supernova remnant", despite 
the shell or blast- wave that would normally correspond 
to a SNR remain i ng as yet unidentifi e d in this source 
( Frail et alJll995t ISeward et all l200d iTziamtzis et al. 
2009). Some PWNe have for historical reasons been 



catal oged as "fiiled-center supernova remnants" (Green) 
2009). while others, discovered more recently, have been 
excluded from such lists. As we describe below, a young 
SNR may contain a PWN, a combination often called 
a composite remnant (e.g., G21.5-0.9; Fig. fT5|) . 



Fig. 14 Images of the Crab Nebula. 



The key differences between SNRs and PWNe are 
as follows: 

— Energy Source: SNRs result from an essentially 
instantaneous deposition of energy, in the form of a 
blast wave driven into the ISM by a supernova ex- 
plosion. In contrast, PWNe have a continuous power 
source, the bulk relativistic flow of electron/positron 
pairs from an energetic neutron star. 

— Radio Morphology: As a direct result of their 
differing sources of energy, SNRs and PWNe have 
distinct radio morphologies. SNRs are usually limb- 
brightened shells of synchrotron emission, while PWNe 
are typically amorphous or filled-center synchrotron 
nebulae, brightest at the pulsar's position. 

— Radio Spectral Index: SNRs usually have rela- 
tively steep radio spectral indices, aw 0.3 — 0.8, as 
shown in Figured! In contrast, PWNe have spectral 
indices in the range, a ~ — 0.3, which is too flat to 
be explained b y simple models of diffusive shock ac- 
celeration (e.g..lAtovanll999t[Fleishman &: Bietenholz 
l2007t rTanaka fc Takaharall2010li ~ 

— Angular Extent: SNRs are long-lived objects with 
a wide range of sizes, with angular extents ranging 
from ~ 1' to > 5°. PWNe are usually relatively 
small, with sizes in the range 10" to 30', although 
a few older PWNe may be significantly larger (e.g., 
Vela X; Fig. EE). 

— Fractional Polarization: At radio frequencies near 
1 GHz, SNRs typically have modest amounts of lin- 
ear polarization, at the level of 5%-10%. PWNe usu- 
ally have very well-organized magnetic fields, with 
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Fig. 1 5 Left: Deep Chandra ima ge of the composite SNR G21.5- 
0.9 jMatheson fc Safi-Harrj|2005h . The PWN is the bright central 
condensation. Right: Schematic structure of a composite SNR 
containing a PWN. 



correspondingly higher polarization fractions, in the 
range 30%-50%. 



2.3 PWN Evolution 

Theoretical expectations predict three main phases in 
the evolution of a PW N (jRevnolds fc Chevalier] 1 1984 : 
Chevalierl[l99il2005h . governed in turn by: 



1. the expansion of the PWN; 

2. interaction of the PWN with the surrounding SNR; 

3. the motion of the pulsar powering the PWN. 



Each of these phases is discussed in more detail be- 



low. 



2.3.1 Phase 1: Expansion into Unshocked Ejecta 

At early stages in a PWN's evolution, the pulsar's spin- 
down luminosity remains relatively constant, so that 
there is a steady injection of energy into the PWN. As 
a consequence, the PWN expands supersonically into 
the surrounding low-pressure environment, with the ra- 
dius, R, of the P WN evolving with time, t, as R oc i 6 / 5 
( Chevalierlll977t) . Thus the PWN drives a shock into 
the inner edge of the expanding SN ejecta (see Fig. ITS]) . 

Many pulsars are born with high space velocities, 
up to and beyond 1000 km s -1 . However, because the 
sound speed inside the PWN is much higher than this, 
the PWN stays relatively centered on the pulsar, and 
has a quasi-spherical appearance. 

Inside the PWN, there are two broad zones. Close 
to the pulsar is the unshocked wind zone, in which par- 
ticles and Poynting flux flow freely outward from the 
central star. The pressure in this wind drops with in- 
creasing distance from the pulsar, until we reach a ter- 
mination shock at which the wind is confined by exter- 
nal pressure. At this shock the wind particles are accel- 
erated up to even higher Lorentz factors, and have their 



pitch angles scattered so that they can emit synchrotron 
radiation. Observations of this emission allow us to in- 
fer the composition of the wind at and immediately 
downstream of the shock. Specifically, the parameter a 
corresponds to the ratio of electromagnetic energy to 
particle energ y at the wind shock, and i s usually much 
less t han one ( Kennel fc Coroniti 1984 : Gaensler et al.l 
2002): the shocked wind is particle-dominated (er -C 1), 
in contrast to the unshocked wind which is expected to 
be radiation-dominated (a ^> 1). 

Beyond the termination shock is the second main 
zone, the emitting region of the PWN, in which the 
bulk flow continues to decelerate to match the external 
boundary conditions, and where synchrotron emission 
is produced. 

The morphology of a PWN in this early phase of 
evolution has been spectacularly revealed by X-ray im- 
ages of the Crab Nebula taken by the Chandra X-ra y 
Observatory (|Weisskopf et al.ll2000t [Hester et alj|2002l) . 
These observations show that the Crab Nebula is dom- 
inated by a bright X-ray torus. The inner part of this 
torus is bounded by a bright X-ray ring, thought to cor- 
respond to the wind termination shock. To the south of 
the pulsar, an X-ray jet runs along the axis of the torus, 
apparently originating very close to the central pulsar 
before curving slightly at large distances from the cen- 
ter. A faint counter-jet can also be seen to the north. 

The torus is thought to result from the Poynting 
flux concentrated in equatorial regions of the system 
by the woun d-up magnetic fiel d. However, a key point, 
first noted bv lLvubarskvl (|2002l ). is that the termination 
shock radius should decrease with increasing angle from 
the equator. The axial jet seen for the Crab Nebula is 

therefore not a jet originating from the pulsar itself, but 

is par t of the post-shock wind (e.g., Boeovalov fc KhangoulvanI 
2002h . 

Relativistic magnetohydrodynamic simulations show 
that at mid-latitudes, the external pressure can reverse 
the post-shock particle flow. The wind flows back along 
the surface of the termination shock toward the poles 
( Komissarov fc Lvubarsky|2003t Ibel Zanna et~alll2004 , 
20061) . at which points it is co llimated into jets by hoop 
stress (the "toothpaste effect: ICamus et al.ll2009l ). Sim- 
ulations thus can do a reasonable job of reproducing 
the morphological properties of the Crab Nebula, as a 
result of a latitude-dependent termination-shock radius 
with a flow reversal. 

Recently published optical polarimetry of the Crab 
Nebula using the Hubble Space Telescope has revealed 
spect acular details of this PWN's large-scale magnetic 
field (jHesterl [2008) . A great deal of fibrous structure 
seen in polarization corresponds to the local field di- 
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Fig. 16 2.4 GHz P arkes image of the Vela supernova remnant 
llDuncan et al.lll996h . The w hite cross indicates t he pulsar, and 
the arrow its proper motion toodson et alj |2003); the fact that 
the pulsar is neither at nor moving away from the PWN's center 
indicates that the reverse-shock interaction has take place. 



rection, while the torus of the nebula appears to be 
confined by a poloidal component of the field. 

2.3.2 Phase 2: Interaction with Reverse Shock 

A PWN moves into a new phase of evolution when its 
outer boundary collides with the reverse shock from 
the SNR in which it is embedded. This typically occurs 
after ~ 7000 years for the typical case of a 10 51 erg su- 
pernova ejecting 10 s olar masses into an ambient I SM of 
density 1 atom cm (|Revnolds fc Chevalierlll984l) . The 
collision with the reverse shock initially compresses the 
PWN, and the system may reverberate for a couple of 
cycles before establishing a new equilibrium (jBlondin et~al 
20011 Ivan der Swaluw et al.ll200ll Ibelfand et al.ll2009h 
The compression of the PWN raises the internal 
magnetic field, resulting in substantial synchrotron burn- 
off at high energies. The X-ray extent of a PWN in 
this phase can therefore be rather small. Furthermore, 
if the pulsar has a high space velocity or if the sur- 
rounding ISM is inhomogeneous, then the reverse shock 
does not collide with all parts of the PWN at the same 
time. As a result of this asymmetric collision, the pul- 
sar can end up substantially offset from the center of 
its PWN, and the PWN's radio morph ology can take 



on a chaotic or filamentary appea rance (Blondin et_al 



200ll Ivan der Swaluw etallboO^ (see Fig.QI]). 



2.3.3 Phase 3: Supersonic Motion 

Since a pulsar moves ballistically while the surround- 
ing SNR decelerates in its expansion, any pulsar with a 
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Fig. 17 Left: Hydrodynamic simulation of a pulsar bow shock. 
The pulsar (asterisk) is moving from right to left with Mach num- 
ber 60. Right: Chandra image (blue) and VLA image (red) of 
G359.23-0.8 2 (the "Mouse"), the bow shock associated with PSR 
J1747-2958 llGaensler et alj 12004^1 . The white arrow indicates a 
structure which may be the termination shock. 



significant space velocity will eventually begin to move 
toward the SNR's rim. For a SNR in the Sedov phase, 
the sound-speed in the remnant's interior dro ps with in- 



creasi ng distance from the center. Specifically. Ivan der Swaluw et al 

(1998) has shown that for a Sedov SNR, the pulsar's 

motion will become supersonic when it has moved 68% 

of the way to the SNR's rim, which occurs at 50% 

of the total time taken for the pulsar to escape the 

SNR completely. Once the pulsar becomes supersonic, 

a new, bow-shock, PWN is formed, which no longer 

expands and which is externally confined by the ram 

pressure resulting from the pulsar's motion. Eventually 

the pulsar escapes the SNR completely, an d now drives 

a bow shock through the ambient ISM ( Bucciantini 

20021 IChatteriee fc Cordesl 120021: iGaensler et al.l 12004 ; 



Bucciantini et al 



2005) (see Fig.QH). 



2.3.4 Measuring PWN magnetic fields 

We can infer the strength of the magnetic field in- 
side PWNe through various techniques: by assuming 
equipartition between particles and magnetic fields av- 
erage d over an assumed synchrotron emitting volume 
(e.g., Irlester et al.lll995l) . by measuring the radius of 



the termination shock from a high-resolution X-ray im- 
age and combining this with a physical model for the 
flow to estimate the field value at this location (e.g., 
IGaensler et al 1l2002l) , or through observations of inverse 
Compton emission in GeV or TeV gamma-rays (e.g., 
MSH 15-55: Abdo et al. 2010b, Aharonian et al. 2005). 
De Jager & Djannati-Ata'i (2008) review GeV/TeV emis- 
sion from PWNe. These approaches show that PWNe 
can have a wide range of nebular magnetic fields, from 
~5 /iG to >1 mG. The magnetic field strength gener- 
ally depends strongly on E, t and c, as well as on the 
pulsar braking index. 
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A variety of magnetic field orientations have been 
inferred in PWNe from the angles of their radio polar- 
ization. Some PWNe show a broadly toroidal magnetic 
field, some show a radial magnetic field, and in others 



there is a complex or tangled appearance. iKothes et al 



observations and modeling of unique systems such as 
PSR J0737-3039A/B (the "do uble pulsar") may pro- 
vide vital insights in this reg ard ( Mc Laughlin et al . 2004; 
Lvutikov fc Thompson! liboa) . 



(2006) lave suggested that all PWNe have predomi- 
nantly toroidal magnetic fields, but that these different 
polarization patterns result from differences in viewing 
angle. Toroidal polarization vectors, such as seen for 
the PWN G106.6+2.9, res ult when the observ er looks 
along the pulsar's spin axis (jKothes et al Radial 
fields, such as seen in the PWN G21.5-0.9 fReichll2002{) 
correspond to a viewing angle perpendicular to the spin 
axis. And more comp lex magnetic field morp hologies, as 
for the PWN 3C 58 dWilson fc Weilerlll976l) . occur for 
oblique viewing angles. The complex polarization prop- 
erties of systems like 3C 58 may also occur as a result 
of kink instabilities, which can produce sheared loops of 
magnetic field th roughout the nebula (|Begelmanlll998t 



Slane et alJl2004h . 



There are only limited observations of magnetic fields 
in pulsar bow shocks. The theoretical expectation is 
that the system forms a toroidal field near the apex of 
the bow shock, but that the magnetic field is collimated 
in the trailing magnetotail. Radio polarimetry matching 
this p rediction has recently been observed by Ng et al.l 
d2010l) for the bow-shock PWN powered by PSR J1509- 
5850 . However, ot her magnetic field geometries have 
also been seen (e.g.. lYusef-Zadeh fc Gaens ler 20051) . sug- 
gesting that the magnetic field structure of bow-shock 
PWNe may depend on whether the pulsar's s pin axis 



3 Constraints on the local interstellar magnetic 
field from radio emission of SN 1006 

3.1 Introduction 

The radio morphology of supernova remnants (SNRs) 
may be very informative of the conditions of the magne- 
tized environments in which the blast-wave expansion 
occurs and, in particular, on the acceleration processes 
at the shock front which give rise to the energetic elec- 
trons ultimately responsible of the synchrotron emis- 
sion in the radio and (possibly) X-ray band. In this con- 
text, the bilaterally symme tr ic or barr e l-shap ed BSNRs 
Kesteven fc Caswelll l|l987l) . iGaenslerl ()l998l )) are con- 



sidered ideal laboratories, because their morphology is 
definitely the result of the lack of small scale inho- 
mogeneties which may confuse the interpretations. A 
point-like explosion in a uniform magnetized medium 
with constant strength and direction of B should give 
rise to a symmetric BSNR whose bright limbs are lo- 
cated where the magnetic field is parallel or perpen- 
dicular to the shock speed, if the injection efficiency 
is greatest where the shock is quasi-parallel or quasi- 
perpendicular/isotropic, respectively, and if B lies in 
ike plane of the sky, whereas d ifferent configurations 



™ „„„„„„ „ " tbe plane 01 the sky, whereas d merent configurations 

aligns or mis-aligns with its velocity vector (e.g.,| Vigelius et ^ at different aspect angles dFulbright fc Reynolds] 



2007). 



2.4 Outstanding Issues and Future Work 

There are many aspects of PWNe which are still active 
areas of investigation. Two particular areas of focus are 
as follows. 

First, several PWNe have now been observed to 
show substantial time variability in their morphology 
and intensity (e.g.JPavlov et al.|[20cililHester et al.l2002 ; 
DeLanev et al.ll2006l : lDe Luca et al.ll2007t l. These obser- 



vations have provided a rich suite of information, but 
we now need relativistic magnetohydrodynamic simula- 
tions to fully understand the flow conditions and insta- 
bilities that produce these effects ( Komiss arov fc Lyubarskvl L _ . .. ,, 
2004 iDel Zanna et al.ll2004 : iBucciantini et al 120051: iBucciam fcriit o the magnetized environments in which the shock 



(Il99i 

However, in real life, BSNRs ar e often asymmetric, 
and therefore Orlando et al. ( 20071 ) ( hereafter Paper I) 
have r ecently generalized the study o fFulbright fc Revnoldsl 
(1990) to the cases in which the explosions occur in a 
large scale gradient of density or magnetic field, show- 
ing that this model is able to reproduce most of the 
asymmetries observed in real BSNRs. In particular, the 
radio morphology loses one axis of symmetry, and the 
limbs are not equally bright (if the gradient runs across 
the limbs) or they are not opposite and they converge 
on the side in which the density or the magnetic field 
is increasing (if the gradient runs between and parallel 
to the limbs). 

It is clear that morphology of BSNRs is tightly cou- 



2010). 



Second, while we have a multitude of observations 
of pulsar winds after they have been shocked, we have 
very limited data on the composition and transporta- 
tion mechanism in the unshocked wind zone. Targeted 



expands, and it is of particular interest to note here 
the preference of BSNR symmetry axes to be oriented 
parall el to the Galactic plane, as reported by iGaensler 
dl998l) . It seems therefore possible to study their mor- 
phology to derive the geometry of the surrounding mag- 
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netic field, thus shedding more light on the microphysics 
of the particle acceleration processes at the shock front. 

The remnant of SN 1006 seems to be the object in 
which this kind of study may be most fruitful. The uni- 
form environment and the bright limbs visible in most 
of the electromagnetic spectrum make it a real case 
study in the field of part icle acceleration mechan isms 
in strong shocks. Indeed. iRothenflug et al.l ([20041 ) . us- 
ing a simple and powerful geometrical argument applied 
to the XMM-Newton X-ray image of SN 1006, based on 
the ratio between the central and the rim luminosity, 
showed that if the remnant is cylindrically symmetric, 
the bright limbs are likely to be polar caps (instead of 
an equatorial belt) and that, therefore, the magnetic 
field is oriented perpendicular to the bright limbs, in 
the NE-SW direction. This argument seems to break 
the dichotomy between the two competing scenarios of 
the dependence of the electron injection efficiency or 
electron acceleration rate on 0Bn, preferring a situa- 
tion in which the injection is most efficient or accelera- 
tion most rapid when the field is along the shock speed 
(quasi-parallel scenario) over the situation in which the 
field is perpendicular (quasi- perpendicular scenario, see 
iFulbright fc Reynolds! ilQffflh l 1 !. 

O n one hand, t hese findings were somehow expected 
since IVolk et al.l (|2003l ) already pointed out that in 
SN 1006 the injection should be maximum at parallel 
shocks. However, in the light of the uncertainties related 
to the details of the acceleration processes, several au- 
thors still conside red the quasi-perpend i cular scenario 
a viable option: Fulbright fe; Reynolds! ( 199dh argued 
against the quasi-parallel scenario pointing out that 
quasi-parallel models often give rise to unobserved mor- 
phologies in the radio band; lYamazaki et al.l (|2004l) still 
considered both models to explain the obs erved width 
of sele cted filam ents of SN 1006 obse r ved b y lBamba et al 
(j2003l) . whereas Amano fc Hoshino ( 2007 ) developed a 
quasi-perpendicular model which agrees very well with 
the same da ta. Moreover, the s ame s imple geometrical 
argument of IRothenflue et al.l (|2004[ ). if applied to the 
radio image, would be in agreement w ith the equatorial 



Table 1 MHD models of SN 1006 used in this work 



belt (cfr. R 



(200 



tt/3 



0.7 in Sect. 3.2 of Rothc nflug et al 



This discrepancy has never been explained in 
the literature and remain one of the most intriguing 



1 Since X-ray synchrotron brightness depends on both the effi- 
ciency of electron injection into the acceleration process, and on 
the rapidity of acceleration to high energies, studies of the X-ray 
morphology involve a combination of both possible effects, while 
radio studies are insensitive to acceleration-rate issues, because 
acceleration to the GeV energies required for radio emission is al- 
ways extremely rapid compared to evolutionary timescales. Thus 
variations of radio morphology with obliquity point to electron 
injection physics alone. 



Name 


V|B^ 


UNIFORM.B 





GRAD1 


1.4 


GRAD2 


1.5 


GRAD3 


2.0 


GRAD4 


4.7 



" The relative variation of the modulus of the magnetic field over 
a scale of 10 pc 



open issues in the comparison between models and ob- 
servations. 

Recently, Petruk et al. ( 2009! ) have further investi- 
gated this issue, showing that, in the framework of a 
simple model of SN 1006 in terms of a point-like explo- 
sion occurring in a uniform density and uniform mag- 
netic field medium, there is no way to reconcile the 
quasi-parallel scenario and the SN 1006 morphology as 
observed in the radio band. This contradiction between 
the radio morphology (suggesting quasi-perpendicular 
injection efficiency scenario) and X-ray morphology (sug- 
gesting quasi-parallel scenario) has prompted us to in- 
vestigate the effects of non-uniformity of the environ- 
ment on the observed properties of th i s remn ant, capi- 
talizing on the work of lOrlando et al.l (|2007l ) on asym- 
metric BSNRs. 

In Sect. 3.2, we will briefly describe the MHD model 
we have used to reproduce the remnant of SN 1006, 
which include a small gradient of the magnetic field. In 
Sect. 3.3, we introduce the methodology for the com- 
parison between modeled and observed radio images of 
SN 1006, while in Sect. 3.4 we will discuss the results of 
the comparison, showing that, not only can the model 
invoking a gradient of |B| reconcile, for the first time, 
the radio and X-ray morphology of the remnant, but it 
also provides stringent constraints on the overall geom- 
etry of the field in the vicinity of the remnant. 



3.2 The model 

Since there is accumulated evidence that the density 
around SN 1006 is fairly constant, with the exception 
of the NW se ctor where an encou n ter with a dense clou d 
is occurring (jAcero et al.l (|2007f ). iMiceli et all |2009t) h 
we argue that a model with a |B| gradient is more ap- 
propriate to describe the asymmetries of the limbs of 
SN 1006. The remnant was modeled as a point-like ex- 
plosion of 1.4 solar masses of ejecta having a kinetic 
energy of E = 1.3 x 10 51 erg occurring in a uniform 
density medium with no = 0.05 cm~ 3 . The ISMF has a 
value of 30 fiG in the environment of the explosion site 
and it is assumed to have a gradient along the Y axis. 
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Fig. 18 Left panel: Definition of 0b ■ Right panel: Definition of 
0V|B|- 

Such a high value (for the location of SN1006 
in the Galaxy) has been chosen in order to take 
into account as a first approximation the effects 
of magnetic-field amplification, to have a post- 
shock |B| of the order of ~ 100/iG in our simu- 
lation, in agreement with observations of some 
SN1006 filaments. Moreover, the magnetic field 
has a gradient along the Y axis. The direction of the 
magnetic field is in the XY plane. V|B| values used in 
our simulations are reported in Table [TJ The gradient in 
the initial conditions has been modeled with a dipole 



locate d along the Y axis, exactly as in I Orlando et al 



(|2007l) . to which the reader is referred for further de- 
tails on the model and numerical code used for the sim- 
ulations. The simulations were stopped at t = 1002 yr, 
checking that the shock velocity and remnant radius 
are compatible with the observed values (4600 km s _1 
and 9.6 pc, respectively). 



3.3 Comparison between models and observations 

3.3.1 The angles defining the viewing geometry 

Several contributions in the literature show that the 
observed morphology of a SNR emitting synchrotron 
radiation is strongly affected by the angle between the 
dominant direction of B and the line of sight, 0b- More- 
over, in our circumstance, since we also have a gradi- 
ent of the magnetic field, we are forced to consider the 
dependence of the observed morphology on the angle 
between the direction of the gradient and the plane of 
the sky 0v|B|- The definitions of 0b and 0v|b| are 
sketched in Figure [T5] For the purposes of this work, 
we have chosen to synthesize our maps in the following 
way: first, we apply a rotation of 0v|B| degrees around 
an axis passing through the center of the remnant and 
parallel to the limbs, where positive angles mean that 
regions of increasing B are closer to us. Next, we ap- 
ply a rotation of 0b degrees around an axis passing 
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O.O0O4 



!l.(*XJ5 



Fig. 19 Radio map of SN 1006 at 1.4 GHz, with the four re- 
gions in which the surface brightness of the remnant has been 
averaged for the purpose of computation of the four morpholog- 
ical parameters: A = RI M max / RI M m i n , C = center / RI M max , 
Rmax = RI M max i / RI M max 2 and O p, the aperture ang le be- 
tween the limb maxima. Adapted from lPetruk et al.l d2009f) . 



through the remnant center and lying in the plane of 
the sky. Others rotation schemes give similar results. 
The adopted values of 0b are from 0° to 90° in steps of 
2° and the ones of 0v|b| are from 0° to 90° in steps of 
15°. Therefore, for each model, we have generated 315 
maps encompassing all the combinations of the relevant 
angles. 



3.3.2 The morphological parameters 

The exploration of the parameter space is a challeng- 
ing task, because it involves the morphological compar- 
ison between many model maps and the real SN 1006 
images. We devised a simple procedure which involves 
the comparison of the value of 4 morphological param- 
eters calculated both from the synthesized radio emis- 
sion maps and from the observed radio image of SN 
1006. The parameters are the ratio between the maxi- 
mum and the minimum surface brightness (Sb) around 
the rim (A), the ratio between the maximum around 
the rim and the center of the remnant (C), the ratio 
of Sb of the two bright limbs {Rmax), and the angular 
separation between the limbs (9 p). The param e ters A , 
Rmax and 6d were also used in lOrlando et al.l ( 2007 ). 
and we refer to that paper for further discussion. The C 
parameter is introduced here to measure the luminos- 
ity contrast between the brightest rim and the center. 
Note that the C param eter does not correspond to the 



Rtt/3 parameter used bv lRothenflug et al.l (|2004l ) to ex- 



clude the equatorial belt scenario for SN 1006, because 
C is measured in a small circular region in the radio 



16 



map (see below). The parameters A, C and R ma x ar e 
measured using the radio map of Petruk et al. (2009), 
by averaging the Sb value in circular regions of 45" ra- 
dius (Fig. |T9j). The values we obtained for the observed 
image of SN 1006 are A = 2.7 ± 0.2, C = 0.36 ± 0.03, 
R max = 1.2 ± 0.1, and 9 D = 135° ± 10°. For the mea- 
surement of the parameters in the synthesized model 
radio maps, we used an automated procedure to find 
the maximum of the two limbs, the minimum between 
rims along the rim and the central position. Then we 
used an average in a circular region whose radius is the 
same percentage of the SNR radius as used in the real 
radio image (5%). Since the model images have been 
synthesized using the same number of pixels per radius 
of the remnant, this procedure ensures that the model 
and observed values of the parameters are comparable. 



3.4 Results 

3.4-1 Quasi-parallel scenario 

In Fig. [201 we show the A - C, R max - C and D - C 
scatter plots computed in the synthesized radio map 
of our SN 1006 model including a weak gradient of 
the magnetic field (B varies by a factor of 1.4 over 
10 pc), and considering a quasi-parallel scenario for 
the obliquity dependence of the electron injection ef- 
ficiency. We also overplotted the values of the param- 
et ers derived us i ng the real radio map of the remnant 



of Petruk et al 



(2009). Though the agreement is not 
exactly perfect in all the plots, we can note that we 
can define a very limited region of the parameter space 
0b — </>v|B| which is compatible with the observed val- 
ues of A, C, Rmax and 9fj. This means that the ob- 
served radio morphology of SN 1006 is overall compat- 
ible with a quasi-parallel scenario for this remnant, if 
we include a weak gradient of B. A nonzero gradient is 
an essential ingredient to reconcile the radio morphol- 
ogy with the quasi-parallel scenario, and this is shown 
in Fig. [21] in which we show the A — C scatter plot 
for various </>b angles, in the uniform B case and quasi- 
parallel, quasi-perpendicular and isotropic scenario^]. 
In this case, this plot suggests that quasi-parallel mod- 
els do not reproduce the observed parameters at all, 
unlike the quasi-perpendi cular and isot r opic s cenarios, 



as already pointed out bv lPetruk et al.1 ([2009). 



Let us now come back to Fig. [20l A remarkable re- 
sult is that the comparison between the quasi-parallel 
model and the observation strongly excludes a situa- 
tion in which the polar caps are in the plane of the 
sky (</>b = 90°) or along the line of sight (0b =0°). 



E 



100 




0.001 



0.010 0.100 

C (center / rim max) 



1.000 



2.5 



2.0 



£ 1.5 



1.0 




0.001 



0.010 0.100 

C (center / rim max) 



1.000 




0.001 0.010 0.100 

C (center / rim max) 



1.000 



Fig. 20 Top panel: A — C scatter plot for our models, assum- 
ing the quasi-parallel scenario and an age of 1000 yr. Different 
colors correspond to different values of </>v|B|> where each dot 
corresponds to a given value of 0b (from 0° to 90° in steps of 
2°. We overplotted the (A,C) values measured for SN 1006 (the 
inset shows a zoom near this point). Middle panel: Same as top 
panel but for the Rmax — O pair. Lower panel: Same as top panel 
but for the 0£> — C pair. 
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Fig. 21 Same as Fig. 1201 top panel, but for the model with 
uniform magnetic field. 



2 In the uniform B case, R„ 



1 and 9 D = 180° always. 
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Fig. 22 Left panel: Synthesized radio image at 1 GHz for our best-fit model (GRAD1 in Table[T] and quasi-parallel injection efficiency 
scenario, </>b = 38° ± 4° and </>v|B| = 15° i 15°), smoothed with a sigma of 2'. Central panel: radio map of SN 1006 of Fig. 1191 but 
smoothed with the same sigma used for the model image. Right panel: same as the left panel, but for the best-fit quasi-perpendicular 
scenario (<^>b ~ 38°). 



The latter geometry would cause a centrally brightened 
m orphology inste a d of t wo limbs, as already discussed 
by lOrlando et al.l (|2007l ) . The best-fit values of the as- 
pect angle we derive from Fig. [50] and a conservative 
estimate of their uncertainties are (ps = 38° ± 4° and 
0v| B | =15° ±15°. 

The synthesized radio map of the best fit model is 
shown in Fig. [22J along with the observed radio map. 
Both maps have been smoothed to hide the small scale 
features and to focus the comparison on large spatial 
scales. In fact, we don't expect a pixel to pixel match 
between model and observations, given the simple as- 
sumptions behind the model. Indeed, the large struc- 
tures of the observed radio emission are very well re- 
covered by the best-fit model and the two images look 
similar. 

The comparison with other gradients, namely the 
GRAD2, GRAD3 and GRAD4 models, shows that it is 
not possible to find a satisfactory fit for all the param- 
eters. In particular, for higher values of the gradient, 
the angle 8d is underestimated, so GRAD1 is the only 
model which gives us an overall good fit. 

3.4-2 Quasi-perpendicular scenario 

We also produced model images in the quasi-perpendicular 
injection scenario, and we repeated the analysis de- 
scribed in the previous paragraphs. In this case, the 
direction of the V|B| is aligned with the direction of B, 
so the angle </>v|b is the same as (/)&, and we did not 
consider it any further. Another consequence of this is 
that Rmax is always 1, so this parameter cannot give 
any diagnostic. The results are shown in Fig. [331 and 
it seems that a good fit can be found for </>b ~ 70° 



(a va lue in agreement wit h the azimuthal profile analy- 
sis of lPetruk etaTl \200$) ). even if the model points in 
the 0d — C diagram seems to be more distant from the 
observed points than any scatter plot in Fig. [201 thus 
indicating a better fit in the quasi-parallel case than the 
quasi-perpendicular case. 



3.5 Discussion 

The proposed method of comparison between models of 
synchrotron radio emission and the real observational 
data of bilateral SNRs is based on the calculation of 
4 morphological parameters, and we have seen that we 
get a good fit for a quasi-parallel injection efficiency 
scenario if we include a gradient of B and a worse fit 
in case of quasi-perpendicular scenario with the same 
gradient. This means that an agreement is found be- 
tween the overall observed morphology of SN 1006 and 
our model, which corresponds to well defined values 
of the viewing geometry angles (0b = 38° ± 4° and 
0v|b| = 15° ± 15°). We know that the angle between 
the SN 1006 axis of symmetry and the Galactic plane 
is roughly 90°, and that its distance is 2.2 kpc and its 
Galactic latitude is 14.6°. Combining all this informa- 
tion, we can plot the direction of the magnetic field and 
its gradient in a 3D representation of the Galactic disk. 
This is shown in Fig. IIMl 

Remarkably, the direction of V|B| points down to- 
ward the plane and the direction of B is aligned with 
the direction of the spiral arm near the remnant. This is 
indeed a very reasonable configuration for the magnetic 
field at this position in the Galaxy, because models of 
Galactic B indicates that at this location above the 
plane the field still retaines spiral arms features, even 
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to be favored by rotation measures of pulsars in the 
plane. At a distance of ~ 550 pc from the plane for 
the SN 1006 remnant, the large-scale field is expected 
to vary by a factor of 1.4 on scales of more than 100 
pc, much longer than the scale-length required by our 
model fit. This opens up the possibility that we are ac- 
tually sampling the rand om magnetic field component. 
Minter fc Spangler (1996) reports length scales for this 



component of the order of few pc, which is compatible 
with the variations we derive. 




0.01 0.10 

C (center / rim max) 

Fig. 23 Top panel: Same as Fig. I20I (top panel) but for quasi- 
perpendicular scenario. There is no <^v|B| angle in this case, be- 
cause B is always aligned with V|B|. Lower panel: Scatter plot 
of 6r> — C parameters. We overplotted the values observed in SN 
1006. 




Fig. 24 Artist's impression of the magnetic field at the location 
of SN 1006 in our Galaxy. The red and green arrows mark the 
directions of the field and its gradient as derived by the best-fit 
model for the synchrotron radio emission of SN 1006, assuming 
a quasi-parallel scenario, derived in this work. The figure is to 
scale. 



if having smaller amplitudes. Therefore, it is tempting 
to conclude that we are sampling the large scale field of 
the Galaxy. However, we note that the best-fit value of 
V|B| (i.e. a variation of 1.4 times over 10 pc) seems to 
be high to be associated with the large scale B of the 
Galaxy. We have used the model of Sun et al. ( 20081) to 
compute the expected variation of the azimuthal field 
in the Galaxy at the location of SN 1006. In particu- 
lar, we have used their ASS+RING model which seems 



3.6 Summary and conclusion 

The synchrotron radio emission of the archetypical bi- 
lateral supernova remnant SN 1006 is compared against 
an MHD and emission model for this remnant, includ- 
ing a gradient of the external magnetic field, parti- 
cle acceleration and its obliquity dependence, namely 
quasi-perpendicular, quasi-parallel, and isotropic sce- 
narios for the injection efficiency. In order to explore 
the parameter space, which is very large due to the ne- 
cessity of including the viewing geometry in the model- 
data comparison, we developed a simplified procedure 
based on the computation of four morphological pa- 
rameters. We have found a very good fit with a model 
assuming quasi-parallel injection efficiency and an as- 
pect angle (/>b = 38° ± 4° between the direction of B 
the line of sight, and </>v|B| = 15° ± 15° between the 
plane of the sky and the direction of the gradient of 
the magnetic field, and a variation of B of about 1.4 
times over a scale of 10 pc. A worse fit is obtained with 
quasi-perpendicular scenario. The overall morphology 
of the observed radio emission at 1.4 GHz is correctly 
recovered by our best-fit model. The projected direction 
of B and V|B| in the Galaxy are along the spiral arm 
and toward the plane respectively, which is in very good 
agreement with the expected direction of the large scale 
Galactic B. However, the implied gradient is too high 
to be associated with the large-scale Galactic B and 
more typical of the random magnetic field components. 

The application of our method to selected samples 
of bilateral supernova remnants may yield independent 
estimates of the geometry of the Galactic field at several 
locations, which can be useful to understand the field 
topology in our Galaxy. 
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